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ABSTRACT 

Reaction rate data have been collected for chemical 
reactions occuring at high temperatures during reentry of 
space vehicles. Large differences in reaction rate data are 
reported by various authors in the literature. 

The principle of detailed balancing is used in modeling 
kinetics of chemical reactions at high temperatures. 
Although this principle does not hold for certain transient 
or incubation times in the initial phase of the reaction, 
it does seem to be valid for the rates of internal energy 
transitions that occur within molecules and atoms. That 
is, for every rate of transition within the internal energy 
states of atoms or molecules, there is an inverse rate that 
is related through an equilibrium expression involving the 
energy difference of the transition. 

Future goals include theoretical modeling to include 
internal enerqy transitions into overall rate expressions 
that can then be used to predict radiation heating effects 
for reentry vehicles under non— equi 1 ibrium conditions. 

NASA Colleague: Carl D. Scott, Ph.D., ED3, X4306 
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INTRODUCTION 


Data from Space Shuttle flights and experimental work 
indicate that non— equi 1 i br i urn radiation due to chemical 
non— equi 1 i br i um in the shock layer of rentry vehicles may 
result in a significant added heat load. (1,2,3) For OTV 
and AQTV with higher reentry velocities, these non- 
equilibrium effects may be even more si gni f i cant . (4 ) 

This work is a continuation of work from last summer which 
will focus in on reaction rate data reported in the 
literature and also consider furthur the principle of 
detailed balancing as applied to transitions within the 
internal energy states of molecules. 

CHEMICAL REACTION RATES 

At an altitude of 80-100 kilometers, the reactant 
substances in the atmosphere include N , 0 , N, 0, NO, 
various positive ions and electrons. However, with these 
few reactants, a multitude of reactions with various 
interdependent relationships could possibly occur upon the 
initiation of a shock wave by a reentry vehicle. The low 
pressure at these altitudes and the high energy input due 
to the shock wave from Mach 20-35 results in a region of 
chemical non-equilibrium between the reentry vehicle's 
surface and the shock front. This chemical non- 


16-4 


equilibrium is not just due to the progress of reactions 
from reactants to products in the ordinary sense of 
chemical reactions. A more significant and complicating 
factor is that internal energies of atoms and molecules 
also take a finite time to change with the rapid input of 
energy. During this time of internal energy equilibration, 
the ordinary rate constant of a chemical reaction is not 
definable. This time of internal energy equi 1 ibration has 
been referred to in the literature as an incubation time or 
a pre-quasi stat i onary time. (5,6,7) 

The chemical reactions that could occur are many and some 
of the representative ones are given on the next page. 

Carlson and Rieper<8) considered only the first twelve 
reactions given on the next page in their investigations on 
medium shock waves (5-7 km/sec). Their numerical 
predictions agreed with their experimental data on the 
electron temperature variation in these shock waves. 
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Chul Park (9, 10) has used all of the above reactions in a 
computer program (NEQAIR) to calculate the extent of 
reactions and radiation in the non-equilibrium zone for a 
shock wave speed of lOkm/sec. At this higher shock speed, 
calculated results did not fit the limited experimental 
data from shock tube experiments. 

There are still di screpanci es in the literature for 
reaction race data as reported by various authors. Some of 
these are given in Table 1 below. 

Some of these differences in reaction rate data are small 
in magnitude and would not affect the results of 
computations drastically. However, many of the reactions 
have orders of magnitude differences in the rate contants 
and also differences in the temperature dependency of the 
reac ^ 1Qns - These include the nitrogen, oxygen and nitric 
oxide dissociations as well as the impact ionization 
reactions. 

Further work on the effect of varying the rate constants 
upon the output of computer programs would seem to be a 
worthwhile effort. Better agreement with experimental data 
may be obtained by adjusting rate parameters within the 
computer program itself. 
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T'abl e 
Rate Constant 
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DETAILED BALANCING PRINCIPLE 


For ordinary, one-step chemical reactions, the principle o-f 
detailed balancing can be shown to -follow -from equilibrium 
relationships. Consider the reaction 
AB + M A + B + M 

at equilibrium, k^ (AB) (M) = k^. (A) (B> (M> 

and therefore k j /k y = <A><B)/<AB> = K, the equilibrium 

constant. 

Even away from equilibrium, k^ /k y = K because the k's 
do not depend upon concentrat i ons. 

Usually kj and kj.are not both available from experimental 
data. Especially at high temperatures, k h is expressed as 
kj /K in computer programs and calculations involving these 
rates of reactions. 

Work last summer led to the conclusion that this 

re i a t i onshi p , kj /k h = K, also called the principle of 
detailed balancing, does not hold under certain conditions. 
If a large, sudden change in temperature occurs in a system, 
such as that generated by a high energy shock wave, then 
the equilibration of energy into the internal states will 
take a finite time and this time may become comparable to 
the reaction time of the system. When this occurs, the 
ordinary chemical rate constant is not well defined since 
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I 
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it usually represents an average rate of reaction over many 
individual rates from internal energy states. There is 
then an incubation time before the reaction rate reaches a 
steady-state or defined value. During this incubation 
time, the principle of detailed balancing does not hold for 
overall chemical reaction rates. (7) 

During the incubation time, internal energy transitions 
became important and the approach to considering all of 
these transitions is summarized in the '' Master Equation 
Approach" which is shown di agramatical ly below. 


MASTER EQUATION APPROACH 



16-10 



This equation could be solved if all k; > , k;; , k/ and k; are 

known. Appro;; i mat i ons to this equation have been tried, 
using detailed balancing for the internal energy 
transitionsik^ and k j- ), and also assuming a harmonic 
oscillator. This development is shown in several references 
such as "Introduction to Physical Gas Dynamics" by Vincenti 
and Kruger. (18) This particular development is for 
vibration levels only and is breifly. summarized below. 
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Probably the best assumption in this development is the 
principle ot detailed balancing -for the internal energy 
states. At the molecular level, where elementary 
collisions cause internal energy transitions, the exchange 
•from tr ansi at 1 onal to vibrational energy or vice-versa in a 
collisional process should be reversible in time by energy 
conservation principles. Thus there is no reason to 
suspect the failing of the principle of detailed balancing 
at this elementary level unless some other energy 
equilibration procedure is occuring that' is not known. 
Rotational energies are lower than vibrational energies, 
but those transitions are assumed to occur very quickly 
even at low pressures so that rotational energies and 
transl ati onal energies are very quickly equi 1 i brated . The 
effect of anharmoni ci ty of the vibrational mode within the 
molecule is an effect that should be considered. This 
effect becomes more pronounced at higher energies and it 
is these higher energies that are important for 
di ssoci ati on . 

Thus the development of the Master Equation using detailed 
balancing and harmonic oscillator assumptions is a valid 
first approximation to calculation of rates of reaction in 
high energy shock waves. However, other approaches are also 
presently being worked on and these are summarised below. 
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OTHER APPROACHES TO NON-EQUILIBRIUM REACTION RATES 


Dove, et al have calculated the effects of coupling for the 
internal modes upon the relaxation of internal energies for 
hydrogen .( 19) They concluded that at high rotational 
energies <J 10), coupling between transl at i onal , 
vibrational and rotational becomes almost as significant as 
translational -roatational coupling and is greater than 
translational-vibrational coupling. Thus at high rotaional 
energies, equilibration in the internal energies occurs 
through combined rotational —vibrational energy level 
changes. 

Koshi has calculated the effect of vi brat i onal —di ssoci ati on 
and dissociation-vibrational coupl i ng (CVDV effect) and 
compared it with experimental data on NO di ssoci at i on . ( 15) 
His agreement is very good and indicates that as the 
translational temperature increases, the vibrational 
temperature lags behind and seems to approach a limit in 
the range of 5000K even though the translational 
temperature is above 7500K. The CVDV effect is due to the 
reduction of the dissociation rate due to the lack of 
vibrational equi 1 i brati on within the molecule and also the 
coupled effect of a slower vibrational equilibration time 
due to dissociation removing energy from the upper 
vibrational levels and thus depleting them too rapidly. 
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This could result in a "bottle-neck" effect. That is, with 
a very rapid and high energy input, such as that from a 
high energy shock wave, the dissociation rate is limited by 
internal energy transitions within vibrational levels and 
not by the dissociation of the molecule from highly excited 
vibrational energies into two atoms. 

Carlson compared experimental results with calculations on 
nitrogen in moderate shock waves studies and included the 
CVDV efeect as well as el ectr on i c-vi br at i onal coupling. 
His claculated electronic temperature and the measured 
values agreed well within experimental error for this 
series of exper i ments. (8) 

Finally, Chul Park has worked extensively on this problem 
and has written a computer program for predicting the non- 
equilibrium radiation in a shock front using available rate 
data. (9,10) His approach utilizes a two or three 
temperature model assuming different temperatures for 
different internal modes of energy and also a translational 
temperature. For some reactions, he uses the translational 
temperature, for some, the vibrational temperature and for 
others, the electronic temperature. Although this is an 
interesting approach, the calculated distances in the shock 
wave for peak and equilibrium temperatues, do not agree 
with experimental data. 
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SUMMARY 


The most significant finding this summer was that there are 
large differences in important reaction rates as reported 
in the literature by various authors. These large 
differences could have significant effects on any computer 
programs or calculations of chemically reacting systems 
utilized in shock wave heat load predictions. It would be 
worthwhile to try some of these other rates in present 
computer programs and check the resulsts with known data. 
The principle of detailed balancing is used widely in 
calculations of reaction rates. It is not valid during 
short incubation times in high energy shock waves in terms 
of the overall dissociation or recombination rate 
expressions. It does appear to be valid at the micro 

energy level concerning transitions within molecules. These 
transitions occur when transl at i onal energy is exchanged 
for vibrational energy in a collision for example. 
However, another significant find this summer is that in 
collisional processes, multi-jump transitions in 
vibrational levels can and do occur. This does not violate 
the law of detailed balancing, but makes it more 

complicated to take this into account. 

Further refinements in calulation on these high energy 
reactions would benefit from the use of coupling effects 
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such as some o-f those mentioned in this paper. These 
include translation-vibration-rotational coupling, 
electronic-vibrational coupling and vibration-dissociation, 
dissociation-vibration coupl ing. 

Future investigations should include models with limiting 
expressions such as the vibrational temperature reaching a 
limit and thus simpli-fy the overall rate expressions. This 
may or may not occur within a time period that is 
insignificant in high energy shock waves. Approaches 
utilizing partition functions and activated complex theory 
also should be investigated. One -final approach -for future 
work would include the information theory approach by such 
authors as Levine to predict products of chemical reactions 
on the basis of excess entropy. (20) 
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